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We propose a method of generating skyrmion vortices in 
a pair of Bose-Einstein condensates occupying two internal 
states of the same atom. We show that a variety of different 
periodic arrays of vortices may be prepared with an appro- 
priate superposition of orthogonal standing electromagnetic 
waves inducing a coherent coupling between the two conden- 
sates. 
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Since the first observations of Bose-Einstein conden- 
sation in dilute atomic gases one important goal 
has been to develop a technology to investigate vortices 
in atomic Bose-Einstein condensates (BECs) . The quan- 
tized vortices in atomic BECs are analogous to the quan- 
tized circulation in superfluid liquid helium Q. In this 
paper we propose a method of generating different peri- 
odic patterns of skyrmion vortices and solitary waves in 
a pair of atomic BECs. The two condensates are coupled 
by electromagnetic (em) transitions. The em fields form 
a standing wave configuration with topological singular- 
ities at zero-field points. The em field amplitude couples 
to the relative phase between the two BECs through the 
atomic transitions. As a result the Rabi oscillations of 
atoms between the two internal states generate an array 
of topological phases in the two condensate system. 

A BEC is expected to exhibit vorticity in a rotating 
harmonic trap if the trap is anisotropic in the plane of 
rotation |^-^. It has also been proposed that stirring 
a BEC in a nonrotating configuration by a laser beam 
p0|-p2| or by an optically-induced potential may gen- 
erate vortices. Topological defects may also emerge as a 
result of a rapid condensation [p^ or in self-interference 



measurements |15 . Recently, it has been suggested that 
vorticity could be imprinted by imaging a BEC through 
an absorption plate ||l^,|l^. Phase- imprinted dark soli- 
tons have been observed in an atomic BEC 

While the previously explained schemes typically in- 
volve a BEC in a single internal atomic state several pa- 
pers [p9|-pT| have also considered the possibility of trans- 
ferring an atomic population from a nonrotating ground 
state to a vortex state by means of a Laguerre-Gaussian 
laser beam consisting of photons with a nonvanishing or- 
bital angular momentum. In a very recent work Williams 
and Holland have proposed a scheme where two in- 
ternal states are coupled and the trapping potentials of 
these states are mechanically rotated. This technique has 
been experimentally realized to demonstrate the forma- 
tion of a vortex in a two-component BEC [p3|. 



In this paper we propose a method of generating vor- 
tices which is related to the previously described schemes 
p9|-p2[ involving internal atomic transitions. We con- 
sider two internal states of the same atom which are 
coupled by em fields via one or multiphoton transitions. 
The em fields consist of two or more orthogonal stand- 
ing waves. With an appropriate wave configuration we 
obtain a spatially-dependent transition strength with an- 
gular momentum singularities at zero-field points. As a 
result of the spatially-dependent driving the internal and 
external dynamics of the atoms are coupled: Around ev- 
ery zero-field point we generate for the two-component 
BEC a topological phase with an integer multiple of 2tt 
phase winding. The proposed scheme has several advan- 
tages: It may be possible to prepare high-quality vortices 
with little noise. The location of an individual vortex 
may be accurately controlled by shifting the zero-field 
point. We may also engineer a large selection of different 
periodic arrays of vortices in the two-condensate system, 
where the periodicity is determined by the driving em 
fields. Moreover, with moving standing em waves we can 
move the positions of the zero-field points, and therefore, 
also the vortex array in a controlled way. 

The dynamics of the two-component BEC occupying 
internal levels |1) and |2) follows from the coupled Gross- 
Pitaevskii equation (CPE) 
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- KiilV^iP + Ki2\'ip2\^) ipi + hn*^2 , (la) 

- 5 + K22|^2P + Kl2|V'in ll>2 + fl^i'l ■ (lb) 



Here the kinetic energy and the trapping potential are 
introduced in H^: 

EE -|1 + ImcoKx' + ajy^ + pjz^) . (2) 

We have also defined the coefficients of the nonlinearities 
iTrti^aijN/m. Here an denotes the intraspecies 
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scattering length in internal level \i) and ai2 stands for 
the interspecies scattering length. The Rabi frequency, 
ri(r), describes the strength of the coupling between the 
two internal levels. The total number of BEC atoms, the 
atomic mass, and the detuning of the em fields from the 
resonance are denoted by N , m, and S, respectively. 

We assume that the Rabi frequency, il(r), is formed by 
two appropriately phase-shifted and orthogonal standing 
em waves, and has the following form: 

f^(i") = -^{sin[fc(x -x)]~i sm[k{y - y)]} , (3) 

where k denotes the wave number of the em field. When 
the phase shifts are equal to zero, x — y — 0, the coupfing 
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vanishes at Xn = Hn = nir/k, for n = 0,±1,±2,- • •. In 
the close neighborhood of the vanishing Rabi frequency, 
\kxn — nirl ^ 1 and {kym — mirl ^ 1, we obtain 

f7(r) (-lr^e-«'^ (4) 



with ^ = (—1)"+'" . Here p = ^ + 'ip- is the radial co- 
ordinate in the xy plane and is the corresponding polar 
angle. In Eq. (|l|) the phase of the Rabi frequency couples 
to the relative phase between the two condensates. The 
phase of ri(r) close to a zero- field point in Eq. (||) has 
the form of the quantized circulation with the unit wind- 
ing number. We show that the coupling between the em 
and matter fields establishes a topological relative phase 
between the two BECs. 

To demonstrate the formation of vortices we numeri- 
cally integrate GPE in two spatial dimensions. For sim- 
plicity, we assume that the traps are isotropic (a^ = 1), 
and that the trapping frequencies and the potential min- 
ima of the two internal states are equal. As an initial 
state we assume a nonrotating ground state for a BEG 
in level |1) and an unoccupied level |2). In the limit of 
strong self-interaction energy the initial state of t/ii may 
be approximated by the Thomas- Fermi solution in 2D: 
V'i(p) = [2(i?2 _ p2)/(7ri?4)]i/2^ for i? > p, and zero oth- 
erwise. Here R = Z[4Kff /(tt^iw/^)]!/* denotes the 2D 
Thomas-Fermi radius of the BEG and I = {■muj)^!'^ is 
the harmonic trap length scale. 

We choose the same ratio between the three scat- 
tering lengths, K^f : : :: 1.03 : 1 : 0.97, 
as for ®^Rb states |F=1, m=-l) and |F=2, m=l) [||. 
In the numerical calculations we use the nonlinearity 
Kl^jihujl'^) = 1000, the Rabi amphtude f^o = 50^/2w, 
and the detuning 5 = 0. As a first example we create a 
single vortex with the unit topological charge or the unit 
quantized circulation. Here the value of the wavelength 
for the em fields A — 201 and x = y = in Eq. (|^), so 
that only one field node is in the condensate. The atom 
population is initially in level |1). The em fields start 
inducing transitions between the two internal levels at 
time t — 0. We decouple the two levels by turning off 
the em fields at time t = 0.02/a; ~ 1.4/r2o. In Fig. |^ we 
display the density \ip2{x,y)\'^ and the phase (l>2 profiles 
of the BEG atoms in level |2) at t = 0.0202/w. We ob- 
serve a vortex in level |2) with a vanishing atom density 
in the center of the vortex core and the 2tt phase winding 
around the vortex line demonstrating the high quality of 
the preparation process. The phase profile of ■(/'i is 
flat corresponding to a nonrotating state. 

The vortex core in Fig. 0is very large due to the atomic 
population in level |1) which occupies the in- 

terior of the vortex and generates a mean-field repulsion 
as also displayed in Fig. Therefore the presence of 
the vortex could possibly be directly verified by imaging 
the density profile of the BEG in level |2) as in Ref. [ p3[ . 
The relative phase between two BEGs has also a dra- 
matic effect on the dynamical structure factor of the two- 



component system p4| , p5[ | which may be observed, e.g., 
via the Bragg spectroscopy p6| , p7|| . 

Next we consider the preparation of several vortices 
by means of the em fields of Eq. (|l|). In Fig. ^ we show 
the density and the phase profiles of rectangular vor- 
tex arrays which are obtained with A = 3/ (left column) 
and X = 51 (right column). The plots are evaluated at 
t = 0.025/a;, when the em fields are still on. The other 
parameters are the same as in the previous example. The 
topological phases are distributed according to Eq. (^) . 

We can investigate the robustness of the prepared vor- 
tex states after the turnoff of the em field coupling. The 
dynamical stability may be studied by integrating GPE 
with 17 = if we ignore the energetic instabilities due to 
the interaction of the BEGs with thermal atoms p8| , ^ . 
In Fig. 1^ we show a vortex at a later time with the same 
set of parameters as in the case of Fig. |[ The em fields 
are turned off at i = 0.02/w. At < = 1.0/w the phase 
has acquired a spiral shape indicating a radial inward 
fiow. Although the BEGs undergo collective radial os- 
cillations, the vortex remains well preserved several trap 
periods. In Ref. the radial oscillations were exper- 
imentally observed. However, after the initial shrinking 
period the vortex expanded and broke up. The 2D coher- 
ent mean-field picture may not be capable of predicting 
the breakup. In the case of a vortex array the dynamical 
deformations are more rapid. However, in our studies we 
have found that the characteristic features of the vortex 
pattern with X = 51 (Fig. ^) still exist even at t = 1.0/w. 

With em fields we can generate a variety of different 
topological structures. If we are able to use transitions 
through intermediate levels involving several photons, we 
may prepare several vortices in arbitrary spatial loca- 
tions with different topological charges. This construc- 
tion would then exhibit a periodicity determined by the 
wavelength of the em fields. In a more general case the 
Rabi frequency from Eq. (^ is determined by 

TL 

m = no n{sin[k, • (r - r(^'))] - ^ sin[q, • (r - r^^'^W^ . 

(5) 

Here we use |kj| — |qj| and kj • q j =0. The exponent pj 
denotes the topological charge of vortex j and n the num- 
ber of vortex arrays. For instance, if we use a two-photon 
transition, we may prepare a vortex with the topological 
charge of two. In this case we only have the j — I term 
withpi = 2 in Eq. (|). 

In Fig. ^ we show a vortex array with the topological 
charge of two at t = 0.02/w, when the em fields are still 
on. In this case A — 51, D,o = SOw, and the other param- 
eters are the same as before. To reduce the phase noise 
in the plottings we set the phase in the figures equal to 
zero when the atom density is very low. Therefore the 
phase in Fig. |^ is zero close to the vortex line. We note 
the increase or the decrease in the value of the phase by 
47r as individual vortex singularities are encircled. 
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In Fig. ^ we have prepared a three- fold symmetric pat- 
tern of three vortices with the unit topological charge. 
This is obtained from Eq. (bh with <\i/q = x, c[2/q = 



0.87y, 



2.7lcij/ 



-0.5i - 0.87j/, qa/g = -0.5£ 

A = 27r/g = 20^, and J7o = 50a;. Butts and Rokhsar 
have shown that a low-energy state of a rotating BEC can 
exhibit a similar symmetry. Nevertheless, the low-energy 
states cannot be directly created even in a rotating trap 
because the rotation has to overcome the energy barrier 
to vortex formation | ■ 

As a final example we prepare dark solitary waves 
with the Rabi frequency ri(r) — ilosin(fca;). In 
Fig.^ wc display the atom density and phase at time 
t ~ 0.025/a;. Here we again use the familiar set of pa- 
rameters, now with \ = bl and f2o = 50[j. 

In the previous examples we only considered stationary 
em fields. We may also obtain a time-dependent Rabi fre- 
quency by generating a moving standing wave from two 
nearly counterpropagating em fields with slightly differ- 
ent frequencies. In that case we can move skyrmion vor- 
tices inside the BECs along with the moving zero-em-field 
points in a controlled way. 

We proposed a method of generating periodic vor- 
tex arrays in a binary BEC. The technique relies on a 
spatially varying standing em wave configuration which 
drives atomic transitions between the BECs. A potential 
experimental limitation is to find em transitions with ap- 
propriate wavelengths to produce vortex arrays of differ- 
ent periodicity and number of vortices. Nevertheless, the 
number of practical transition frequencies may possibly 
be increased also by tuning the strength of the trapping 
potential. 

In the present work we have demonstrated only in 
terms of simple examples the formation of vortices and 
solitary waves. More complicated em field configurations 
could possibly generate more complex topological struc- 
tures. The proposed technique could obviously be also 
generalized to multi-component BECs. 

We acknowledge discussions with J. R. Anglin. This 
work was financially supported by the EC through the 
TMR Network ERBFMRXCT96-0066. 




FIG. 1. The preparation of a single vortex with the unit 
circular quantization. Before the driving em fields are turned 
on, a nonrotational BEC occupies level The em fields 
generate a vortex in level |2). We show the density profile 
\'tlj2{x,y)\^ (upper left) of atoms in level |2), the density along 
the X axis (upper right) in levels |2) (solid line) and |1) (dashed 
line), and the phase profiles of levels |2), (j)2, (lower left) and 
(f)!, (lower right) at a later time. The density and phase 
of level 1 2) remarkably clearly display the characteristic prop- 
erties of a unit-quantized vortex with the 2iv phase winding 
around the vortex core and the vanishing atom density in the 
center of the core. There is no vorticity in level 1 1) . 
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FIG. 2. The preparation of rectangular arrays of vortices 
with the unit circular quantization. We show the density 
\'tp2(x,y)'\^ and the phase 02 profiles of level |2), and the den- 
sity along the x axis for atoms in levels |2) (solid line) and |1) 
(dashed line) for A = 3/ (left column) and for \ = hi (right 
column). The circulation around every individual vortex is 
associated with a 2n change of phase. The flow around the 
neighboring singularities exhibits opposite helicity. 
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FIG. 3. The dynamics of a single vortex. We show the 
phase 02 of atoms in level |2) at t = 1.0/lj with the same 
parameters as in Fig. |l| The spiral shape of the phase indi- 
cates inward flow due to radial oscillations. The atom density 
is displayed in levels |2) (solid line) and |f) (dashed line) at 
t = lA/u, when the BEC in |2) starts expanding, and in levels 
|2) (dashed-dotted line) and jl) (dotted line) aX t = 2.5/a;. 
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FIG. 4. The preparation of a rectangular array of vortices 
with the topological charge of two. We show the density and 
the phase profiles of level |2) for A = 51. The phase profile 
displays the 47r phase windings around the individual vortex 
lines. 
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FIG. 5. The preparation of a three-fold symmetric pattern 
of vortices with the unit circular quantization. We display the 
density and the phase of atoms in level |2). The circulation 
around the whole pattern changes the phase by 6tt. 
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FIG. 6. The preparation of an array of dark solitary waves. 
We show the density and the phase in levels |2) (solid line) and 
|1) (dashed line) along the x axis for A = 51. The phase graph 
displays a sharp phase slip of the order of tt in the center of the 
solitary wave corresponding to the vanishing atom density. 
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